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The MHD Kelvin-Helmholtz Instability II: 

The Roles of Weak and Oblique Fields in Planar Flows 

T. W. Jones Joseph B. Gaalaas^’^’”^, Dongsu Ryu and Adam Frank^’®’® 


ABSTRACT 


We have carried out high resolution MHD simulations of the nonlinear evolution of 
Kelvin-Helmholtz unstable flows in 2^ dimensions. The modeled flows and fields were 
initially uniform except for a thin shear layer with a hyperbolic tangent velocity profile 
and a small, normal mode perturbation. These simulations extend work by [Frank e\ 
al. (1996) and [Malagoli, Bodo &: Rosner (1996). They consider periodic sections of 


flows containing magnetic fields parallel to the shear layer, but projecting over a full 
range of angles with respect to the flow vectors. They are intended as preparation 
for fully 3D calculations and to address two specihc questions raised in earlier work: 
1) What role, if any, does the orientation of the field play in nonlinear evolution of 
the MHD Kelvin-Helmholtz instability in 2^D. 2) Given that the field is too weak to 
stabilize against a linear perturbation of the flow, how does the nonlinear evolution 
of the instability depend on strength of the field. The magnetic field component in 
the third direction contributes only through minor pressure contributions, so the flows 
are essentially 2D. In Frank et al. (199^ we found that fields too weak to stabilize a 
linear perturbation may still be able to alter fundamentally the flow so that it evolves 
from the classical “Gat’s Eye” vortex expected in gasdynamics into a marginally 
stable, broad laminar shear layer. In that process the magnetic held plays the role of a 
catalyst, briehy storing energy and then returning it to the plasma during reconnection 
events that lead to dynamical alignment between magnetic held and how vectors. In 
our new work we identify another transformation in the how evolution for helds below 
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a critical strength. That we found to be ~ 10% of the critical field needed for linear 
stabilization in the cases we studied. In this “very weak field” regime, the role of 
the magnetic field is to enhance the rate of energy dissipation within and around the 
Cat’s Eye vortex, not to disrupt it. The presence of even a very weak field can add 
substantially to the rate at which flow kinetic energy is dissipated. 

In all of the cases we studied magnetic field amplification by stretching in the vortex 
is limited by tearing mode, “fast” reconnection events that isolate and then destroy 
magnetic flux islands within the vortex and relax the fields outside the vortex. If the 
magnetic tension developed prior to reconnection is comparable to Reynolds stresses 
in the flow, that flow is reorganized during reconnection. Otherwise, the primary 
influence on the plasma is generation of entropy. The effective expulsion of flux from 
the vortex is very similar to that shown by [Weiss (1966)| for passive fields in idealized 
vortices with large magnetic Reynolds numbers. We demonstrated that this expulsion 
cannot be interpreted as a direct consequence of steady, resistive diffusion, but must 
be seen as a consequence of unsteady fast reconnection. 


Subject headings: instabilities - magnetic fields - magnetohydrodynamics: MHD - 
turbulence 
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1. Introduction 


Weak magnetic fields threading conducting fluid media can play vital dynamical roles, even 
when traditional criteria, such as relative magnetic and gas pressures, suggest the fields are 
entirely negligible. Perhaps the best example of this is the destabilizing influence of a vanishingly 


small magnetic field crossing a Keplerian accretion disk ( Balbus &: Hawley (1991) ), where the 
mere presence of the field seems fundamentally to alter the local flow properties. Other examples 
abound, however, that could be particularly important in astrophysics. Among them, we would 
include weak fields penetrating turbulent or otherwise strongly unstable flows, where such fields 


significantly alter evolution and transport properties, {e.g., Biskamp fc Welter (1989) ; Cattaneo 
I&: Vainstein (199T] ; Nordlund et al. (1992) ; Jun, Norman &: Stone (1995) ). Sheared motion is a 
critical common element in many of these flows, and the consequent stretching of a weak, but 
large-scale field can lead to a locally enhanced role for the field. The flows will also frequently lead 
to current sheets and associated magnetic field topologies unstable to reconnection, and that is 
central to the nonlinear evolution of the systems {e.g., Biskamp (1993)| ; Parker (1994)| ). Through 
these processes the fields can also have more global consequences. Study of the nonlinear evolution 
of the classical Kelvin-Helmholtz (KH) instability could be particularly useful as a well-defined 
example of strongly sheared flows. Further, since KH unstable boundary layers are probably 
common, the behavior of the instability is important for its own sake. 


Although the KH instability is fairly well studied in ordinary hydrodynamics {e.g., Corcos| 


fc Sherman (1984) ), comparable study has been much slower in magnetohydrodynamics (MHD). 
That is because the magnetic field substantially complicates the physics itself and also because 
computational methods and resources needed for such studies are only recently up to the task. 
The linear analysis of the MHD KH instability is relatively straightforward and was long ago 
carried out for a number of simple flow and field configurations {e.g., Chandrasekhar (1961); 


Miura &: Pritchett (1982) ). Generally, and especially if the velocity change is not supersonic, the 
ordinary fluid shear layer is unstable to perturbations with wave vectors in the plane of the shear 


layer and with wavelengths greater than the thickness of the layer {e.g., Miura (1990) ). When 
there is a field component projecting onto the flow field, magnetic tension provides a stabilizing 
influence. A simple vortex sheet is stabilized against linear perturbations whenever the magnetic 
field strength is sufficient that Cq > |(ic • Uo)/{2k ■ Bo)|, where Uq is the velocity difference between 
the two layers, Ca is the Alfven speed, k is the perturbation wave vector, and Bq is the direction 
of the magnetic field ( [Chandrasekhar (1961) ). 

I I 


Frank et al. (1996) (Paper I) and Malagoli, Bodo &: Rosner (1996)| (MBR) recently presented 
complementary nonlinear analyses of the MHD KH instability in mildly compressible flows based on 
two-dimensional numerical simulations carried out with new (and different) Riemann-solver-based 
MHD codes. While not the first numerical studies of the MHD KH instability, they represented 
big improvements over previous calculations in both numerical resolution and extent to which flow 
evolution was followed towards asymptotic states (readers are referred to Paper I for additional, 
earlier citations). Considering perturbed “2D” flows that were uniform except for a thin, smooth 
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velocity transition layer, those two papers emphasized the qualitatively different behaviors in the 
nonlinear evolution of unstable flows depending on how close the field strength is to its critical 
strength for stabilization. For fields only slightly below the critical value, enhancements in the 
tension of the field through linear growth can stabilize the ffow before it develops distinctly 
nonlinear characteristics. For weaker fields, however, the initial evolution of the instability is very 
similar to that for the ordinary KH instability. That results in the formation of eddies, and hence 
to substantial stretching of the magnetic field lines as well as reconnection. Paper I emphasized 
the remarkable fact that in a case with a field 2.5 times weaker than critical, reconnection can 
lead to self-organization in the flow and fairly rapid relaxation to a quasi-steady laminar and 
marginally stable flow. MBR presented summaries of simulations extending to somewhat weaker 
fields showing evidence for similar behaviors. 

Neither Paper I nor MBR, however, explored the problem in sufficient depth to establish 
the conditions necessary for the previously mentioned self-organization. In addition, it is very 
important in this situation to understand how the magnetic fields behave when they are “very” 
weak (a concept whose definition needs clarification, in fact). A closely related matter is what 
differences, if any, exist between the behavior of a truly weak field and a stronger field whose 
projection onto the ffow vectors is weak. Alternatively stated, are there differences between the 
nonlinear “2D” MHD KH instability and the “2^0” MHD KH instability? Answers to those 
basic questions are the objective of this paper. We find: 1) for the cases we have considered with 
an initially uniform field that the magnetic field transverse to the plane is unimportant and 2) 
there is a transition from the role of the magnetic field as a catalyst to flow self-organization to 
a role as an added source of energy dissipation that should vanish directly as the initial magnetic 
field strength projected onto the plane vanishes. Ultimately we must understand the full “3D” 
problem, in which the perturbation wave vector also lies outside the ffow direction. On the other 
hand it has been difficult to carry out 3D MHD simulations with sufficient numerical resolution to 
be confident of the results in complex flows such as these. In addition, it will be useful to compare 
fully 3D behaviors with 2D flows. We hope that the current work is a significant, constructive 
step towards a full understanding of this problem. 

The paper plan is as follows. In §2 we will summarize the problem set-up and relevant results 
from Paper I. §3 contains a discussion of new results, while §4 provides a brief summary and 
conclusion. We also include an appendix presenting an analytical model for diffusive flux expulsion 
from a steady vortex, in order to contrast that physics with what we observe in the eddies that 
form in our simulations. 


2. Background 

In order to focus on specific, important physical issues we have chosen to explore an 
idealization of the MHD KH problem, reserving for the future the more general problem. We 
present in this section only a bare outline of our method and some key results from previous work. 
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A full discussion of the computational setup along with several tests of such issues as adequate 
numerical resolution and geometry of the computational box can be found in Paper I. 


2.1. Problem Definition 


The geometry of the computations is shown in Figure 1. The only difference from Paper I is 
that in those earlier computations we assumed an aligned field, 9 = cos“^ \Bo ■ Uo\ = 0, whereas we 
now relax that constraint to include magnetic fields oblique or orthogonal to the flow plane. We 
assume the flow to be periodic in x and that the y boundaries are reflecting {i.e., neither flow nor 
field lines cross the y boundaries). This was a configuration used initially by Miura (1984), and 
we followed it in Paper I to enable a direct comparison with his results. The influences of those 
boundary choices are discussed fully in Paper I. In brief, the periodic boundary limits coalescence 
of structures to scales equal to the box dimension, L. That is more significant than the existence 
of the reflecting boundaries, which seem to have only minor influence on dynamics in the (narrow) 
central regions where flow organization and dissipation is largely determined. 


The initial background flow has uniform density, p = 1, gas pressure, p = 0.6, and an 
adiabatic index, 7 = 5/3, so that the sound speed, Cg = \/"ip/p = 1-0. The magnetic field, 

Bq = Bo{x cos 9 + zsin0), is also uniform. In Paper I we considered cases with Bq = 0.4, 0.2 (and 
0 = 0), so that Ma = 2.5, 5, since those were studied by [Miura (1984) . We now add to these a 
number of new cases as outlined in Table 1. To facilitate comparisons we identify the simulations 
from Paper I as Cases 1 and 2 in Table 1, with the new simulations following. The velocity in the 
background state is antisymmetric about y = L/2 according to the relation 


Uo = Uo{y)x = - ^ tanh 

with Uo = 1- This describes a smoothly varying flow within a shear layer of full width 2a. For 
all our simulations presented here a = T/25, chosen to make the interactions with the reflecting 
boundaries small. The square computational box has L = 2.51. Flow is to the left in the 
top half-plane and to the right below that. To this state we add a perturbation, 5{p,p, u,B), 
defined to be a normal mode found from the linearized MHD equations appropriate for the 
chosen background, periodic in x and evanescent in y, with period equal to the length of the 
computational box, L. This was done exactly as in Paper I. We note that when 9^0 there 
are perturbations in all three vector components of the magnetic and velocity fields. Under the 
conditions we used the computational frame is comoving with the KH waves. All flow and field 
quantities are either symmetric or antisymmetric around two points, which happen with our 
choice of phases in 5{p,p, u,B) to be at y = T/2, x = T/4, ‘iUjA. Because the velocity field is 
antisymmetric around these points they are the places where strong vortices tend to form. 


The equations we solve numerically are those of ideal compressible MHD; namely, 

| + v.M = o, 


( 2 - 2 ) 
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d u 11 

^ + u • Vu H — Vp {V X B) X B = 0, 


dt 


P 


P 


dp ^ ^ 

+ u-Vp + 7 pV • u = 0, 

dB 

- V X (u X B) = 0, 


(2-3) 

(2-4) 

(2-5) 


along with the constraint V • B = 0 imposed to account for the absence of magnetic monopoles 
(e.g., Priest (1984) ). The isentropic gas equation of state is p oa p'^. Standard symbols are used 
for common quantities. Here, we have chosen rationalized units for the magnetic field so that the 
magnetic pressure pb = B^/2 and the Alfven speed is simply Ca = B/y^. 

These equations were solved using a multidimensional MHD code based on the explicit, finite 
difference “Total Variation Diminishing” or “TVD” scheme. The method is an MHD extension of 


the second-order finite-difference, upwinded, conservative gasdynamics scheme of Harten (1983) 


as described by Ryu &: Jones (1995). The multidimensional version of the code, along with a 


description of various one and two-dimensional flow tests is contained in Ryu, Jones &: Frank 


(1995) . The code contains an j^t-based routine that maintains the V • B = 0 condition at each time 
step within machine accuracy. That step does not compromise the other conservation relations. 


Numerical solution of equations (2-2) - ( p-5| ) on a discrete grid leads, through truncation 
errors, to diffusion of energy, and momentum, as well as to entropy generation. Of course, 
such effects are also present in nature and are important to defining the character of the flows. 
The existence of effective numerical resistivity is necessary, for example, to allow magnetic 
reconnection to occur in the calculations. Our methods exactly conserve total energy, as well 
as mass, momentum and magnetic flux, so the exchange between kinetic, thermal and magnetic 
energies along with entropy production is internally consistent. There is fairly good evidence 
that conservative monotonic schemes, as this one is, do a good job of approximately representing 
physical viscous and resistive dissipative processes that are expected to take place on scales 


smaller than the grid {e.g., Porter &: Woodward (1994) ). For the astrophysical environments 
being simulated the expected dissipative scales are likely very much smaller than those that can 
be modeled directly. Recent numerical studies of reconnection suggest in the MHD limit with 
large kinetic and magnetic Reynolds numbers that the local energy dissipation rate through 
reconnection becomes independent of the value of the resistivity in complex flows {e.g., [Biskamp 


(1993); Biskamp (1994)1 ). However, when we depend on numerical dissipation, we must be cautious 
about the possible role of “uncaptured” dynamical structures that could be expected on scales 
smaller than the grid scale, or about magnetic field structures on such scales that could enhance 
the number of reconnection sites. In fact, we shall see that for our KH instability-induced flows 
involving weak fields (where reconnection topologies are formed on many scales) total energy 
dissipation is slightly “enhanced” when the grid is finer. This is opposite to what we expect from 
the effects of reduced numerical diffusion alone, but consistent with an increase in the number 
of reconnection sites, or “X” points, allowed when smaller-scale field structures can be resolved 
cleanly.. Further, we find that for a fixed numerical resolution, addition of a very weak magnetic 
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field substantially enhances the rate of energy dissipation, again as one expects in response to 
reconnection {e.g., dimmer, Lesch &: Birk (1997) ). Thus, we see strong evidence for unsteady, 
local reconnection as in high Reynolds number, MHD turbulence, but also that our numerical 
solutions are not quite converged in terms of total dissipation. 

Below we will express all of our results in time units defined by the growth time of the linear 
instability, tg = T”^, as estimated from graphs presented by pVIiura &: Pritchett (1982) . That is, 
we express time as r = t/tg. We find that especially the initial saturation of the instability, but 
also the relaxation processes are fairly uniformly expressed in these units. The time units, tg, are 
listed for each case in Table 1, along with the duration of the simulation in units of r. The flows 
examined in the present paper all have tg ~ 1.6 — 1.7. For comparison the sound crossing time in 
the box IS tg = L = 2.51 (since Cg = 1) and the Alfven wave crossing time is = M^tg (since 
Uo = 1). The normalized turnover time for a large eddy is roughly tE L/{U,/2) 2tg. 


2.2. Magnetic Field Evolution 


Since it enters prominently into our later discussions, it is helpful here to remind readers in a 
simple way of what we can expect for the local evolution of the magnetic field in these simulations. 
A full analysis of all the subtleties is beyond our scope here, so readers are referred to detailed 
discussions such as those by Moffatt (1978) , Priest (1984) , and Biskamp (1993)| . Equation (2-5) 
follows fields when the resistivity is exactly zero. As already mentioned, our finite difference code 
will introduce effects that mimic a hnite resistivity, although it is not possible to define an exact 
value for the resistivity, rj. The effective resistivity will also depend on grid resolution, decreasing 
roughly as N~‘^ within smooth flows. Despite these limitations we can make good heuristic use of 


the resistive MHD extension of the induction equation (^). For our purposes it is interesting to 
cast that equation in the following form: 


d ln(|B|/p) 
dt 


1 d Inph 

2 dt 


d \np B ■ [(B • V)u] — -k 7?V • (j x B) 


dt 


B2 


( 2 - 6 ) 


where d/dt is the Lagrangian time derivative, the first term on the right represents field 
amplification by stretching, and the last two terms containing the resistivity account for magnetic 
“annihilation” and “diffusion”. We have also used equation |2-2| . The current density, j = V x B. 
The term is a dissipative term that balances the Joule heating in the analogous energy 
equation for the gas {e.g., Paper I). The last term is written as an expression involving the Lorentz 
force, j X B, to show that it represents the transport of momentum flux in response to resistivity; 
i.e., the slippage of field lines. 


A frozen-in field results when rj = 0. That leaves only the time derivatives and one term on 
the right of equation [ 


If there is flow compression or expansion only perpendicular to B then 
the right side of equation 2-6 vanishes and leads to \B\/p = constant or (x p^. Those are the 
most common statements of field compression. However, more generally one needs to include the 























other ideal MHD term on the right that accounts for “stretching”. In fact, field enhancements 
due to compression are much more limited than those due to stretching, especially in mildly 
compressive flows, such as those we are studying. 


We also can see from equation (2-6) that resistive influences on magnetic energy are associated 
with both Joule heating and momentum transport. In fact reconnection leads to both irreversible 
heating of the local plasma and to its acceleration. The physics of reconnection is complex and 
beyond the scope of this paper. However, it may be helpful to estimate the dissipation rate using 
equation |2-6| and the Sweet-Parker description of reconnection {e.g., Biskamp (1994 )|; Parker 


(1994) ; Zimmer, Lesch Sz Birk (1997)1) . Assuming an incompressible flow steadily carrying 


oppositely directed fields into a current sheet of thickness, 5, equation ^ gives us a magnetic 
energy annihilation rate and associated dissipation rate per unit volume, Q = ~ where 

u is the inflow speed. In this picture plasma flows out from the reconnection region at the Alfven 
speed, so mass conservation leads to the relation between the current sheet thickness, S, and 
its width, 1 = 5^. Then, Q ~ So, the integrated dissipation rate through reconnection 

depends on the field energy advected into reconnection sites and the summed volumes of all the 
reconnection sites. The current density within the current sheet can be estimated as j ~ B/5, so 
that the aspect ratio of the reconnection region is 5/1 ^ l/y/Ni, where Nl = {I Ca)/r] is known as 
the Lundquist number of the plasma and is obviously related to the magnetic Reynolds number. 
Consequently the dissipative volume scales as t^/y/Wj/. In 2D, dissipative reconnection regions 
form out of tearing mode instabilities within a current sheet when the aspect ratio, 5/1, is small 
( Biskamp (1994) ), and, hence when the magnetic Reynolds number is large. Thus, reconnection is 
not really steady, and the number of sites and their individual volumes will depend on the magnetic 
Reynolds numbers in the flow. Those are, indeed, the behaviors we see in our simulations. 


2.3. Issues 


We already alluded in the introduction to the basic character of nonlinear MHD KH instability 
properties found from previous work. Our intent here is to explore more fully the behaviors of 
weak magnetic helds in this situation. In preparation for that we note from past work several key 
features for 2D symmetry: 


1) When there is no magnetic field or if the field is orthogonal to the flow direction a shear 
layer will “role up”. For the periodic flows considered here the result is a stable “cat’s eye” vortex 
whose length equals the imposed periodicity on the space and whose height ~ 1/3 the length. As 
long as the flow is subsonic or “submagnetosonic” no shocks are involved and the vortex decays 
only through viscous diffusion. We will not deal with the supersonic or supermagnetosonic cases 
here (for some work on those see, e.g., Pedelty &: Woodward (1991) ; Miura (1990) and references 
therein). 


In many astrophysical applications the kinetic Reynolds number of the flow is very large, so 
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we wish to consider similar cases, so that dissipative decay times are long. For our simulations, 
the empirical viscous decay time of the flow is at least four orders of magnitude longer than 
the duration of our computations (see Case 5 energy evolution curves in Figure 3). Thus, the 
cat’s eye represents a “quasi-steady relaxed state”, in which the shear layer has spread vertically 
by horizontal localization of vorticity and become stable. Of course, for non-periodic systems 
there will be continued spreading due to additional vortex mergers, while in 3D the flows will be 
unstable to perturbations directed along the third direction. Those influences are beyond the 
scope of our present investigation, however. 

2) In the other extreme, if the magnetic tension force produced by a perturbation in the shear 
layer exceeds the “lift” force produced by the perturbation, the perturbed flow is stabilized. For 
wavevectors aligned with the flow that condition exists in a linear perturbation whenever the field 
is strong enough that = Uo/{BoCos6) < Mac = 2. Then only viscous diffusion contributes 
to spreading of the original shear layer. As before, we can neglect that influence over finite times 
and describe the shear layer as remaining in a quasi-steady relaxed state from the start. For the 
present computations the critical magnetic field to stabilize a linear perturbation is Be = 0.5. 

For fields slightly weaker than critical, a small, but finite amplitude perturbation may lead 
to the same stabilization condition, perhaps after a small amount of quasi-linear growth to the 
instability. Such was the result obtained for the so-called “strong-held” case of Paper I, where 
Ma = 2.5. That calculation is listed as Case 1 in Table 1, here. Under these conditions the flow is 
never far from laminar and there is a modest amount of spreading in the shear layer before it also 
reaches a quasi-steady relaxed state consisting of a broadened, laminar shear layer (a point made 
both in Paper I and by MBR). In either Case 1 or Case 2, the total magnetic energy changed very 
little during the flow evolution. In Paper I we pointed out that for the symmetry imposed the mean 
vector magnetic held is time invariant; i.e., {B) = constant, so any relaxed state with a relatively 
uniform held will automatically contain magnetic energy close to that of the initial conditions. 
This condition just rehects the conservation of magnetic hux on the grid. The “relaxed” magnetic 
energy is very slightly enhanced, because the held is not quite uniform at the end. A small amount 
of kinetic energy dissipation takes place in Case 1 before the how becomes relaxed. That amount 
is mandated by total energy conservation, the approximate magnetic energy conservation in this 
case, and the fact that on a hxed space the kinetic energy in a broad, symmetric shear layer is 
less than in a thin one. So, any evolution leading to a broadened, laminar shear layer requires an 
increased thermal energy determined by the hnal width of the layer, independent of how it got 
there. In our idealization the mass in the box is also exactly conserved, so the mean density is 
constant. Consequently, to hrst order there is no change in thermal energy by way of reversible, 
adiabatic processes. Most of the increased thermal energy must result from entropy-producing, 
dissipation of some kind. Putting it simply, for these hows to relax entropy must be generated. 
These energy considerations apply to all of the calculations in our study if the relaxed state is a 
laminar how. (It turns out that the kinetic energy of the cat’s eye vortex is less than that for the 
initial how, as well. So, its formation must also generate entropy.) 
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3) For initial fields too weak to prevent formation of the cat’s eye vortex through magnetic 
tension, the shear layer will role up as for unmagnetized fluid flow. In Paper I, however, we saw 
that when the initial field is only a few times weaker than the critical value for linear stabilization, 
there follows a dramatic transformation in the flow as the cat’s eye develops. Paper I considered 
a flow with Ma = 5. We list it in the current Table 1 as Case 2. As the vortex roles up in such 
cases, magnetic field lines are stretched around it, thus increasing magnetic energy at the expense 
of kinetic energy of the flow. The greatest magnetic pressures are produced in thin flux tubes 
formed between the vortex and its twins in periodic extensions of the space (see Figures 4a &: 6). 
Especially within the vortex and around its perimeter, this evolution leads to magnetic reversals, 
beginning after about one turn-over time for the vortex. The reversed helds are unstable to tearing 
mode reconnection. So, once that happens the magnetic field quickly reorganizes itself, releasing 
the magnetic stresses and stored magnetic energy. For Case 2 this leads, as well, to disruption of 
the cat’s eye and eventually to an almost steady, laminar flow, after formation and disruption of 
several weaker vortices. That, final, quasi-steady relaxed state, was similar to the initial conditions, 
except that the shear layer was broad enough to be stable against perturbations on scales that 
fit within the periodic box. The relaxed shear layer had a linear velocity profile (see MBR for 
similar points). There are additional interesting characteristics of the relaxed state. Remaining 
fluctuations in the magnetic and velocity fields were almost exactly correlated, so that they could 
be described as linearly polarized Alfven waves. The magnetic energy returned to a level slightly 
above the initial conditions, with the final excess representing a pair of apparent “magnetic flux 
tubes” bounding a hot central core of the shear layer containing most of the entropy generated 
during the relaxation process. This final condition was reached by about r = 20. The relaxed 
shear layer was broader in this case than Case 1 (see also MBR), so that the kinetic energy was 
also smaller. Thus, even accounting for the slight increase in magnetic energy, we could correctly 
predict that Case 2 with a weaker field was necessarily more dissipative than Case 1. That extra 
dissipation could come about only through the effects of reconnection. Other aspects of this case 
will be visited presently, since we will encounter them again in some new cases examined here. 

Thus, from the calculations reported in Paper I and in MBR it is obvious that weak magnetic 
fields can play a major role in the evolution of MHD KH unstable flows. It is also apparent that 
these roles involve the exchange of energy and momentum from the gas to the magnetic held and 
then back to the gas through Maxwell stresses and reconnection. But, it is not yet clear what 
are the crucial steps in that exchange, nor how it depends on the initial strength of the held, so 
that the ordinary KH behavior results, if it does, in the limit that the magnetic held becomes 
vanishingly small. In addition, although the linear MHD KH instability is not affected by the 
presence of a (possibly strong) held transverse to the how (but still aligned to the plane of the 
shear), a held oblique to the how can carry “circularly” polarized Alfven waves. Since we found 
that the 2D version of the problem generated linearly polarized Alfven waves, it may be important 
to see if the nonlinear problem depends at all on the orientation of the held, or just on the strength 
of the held projected onto the plane as in the linear problem. 
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3. Results 

To meet the objectives at the end of the previous section we have carried out a set of seven 
new simulations. They are outlined in Table 1 as Cases 3-9. (Once again, Cases 1 and 2 were 
discussed in Paper I.) The new simulations were designed to cover a wide range of strengths for 
the magnetic field in the computational plane. They include flows that have exactly the same 
total field strength as Case 2, but in which the initial field is oblique to the computational plane; 

1. e., 0 7 ^ 0 (Cases 3, 5 and 6), as well as flows in which the field is entirely in the computational 
plane, but is weaker than that in Case 2 (namely Cases 4, 7, 8 and 9). The plane-projected field 
strength for Cases 6 and 7 is an order of magnitude weaker than for Case 2. Note that Cases 3 
and 4 are paired to have the same initial planar field strengths, as are Cases 6 and 7. This enables 
us to compare efficiently any distinct roles of field strength and orientation with respect to the 
plane. In Case 5 the field is orthogonal to the plane, so we expect (and see) no important role 
for the magnetic field. As others have noted before, compressible influences are controlled in that 
case by magnetosonic waves rather than pure sound waves, so there is a very slight modification 
in response to that (M“^ = \j -|- M~^ = 1.02 in Case 5). We will not concern ourselves at 
all with flows in which the field is strong, by which we mean situations where magnetic tension 
precludes the nonlinear development of the MHD KH instability. 

Figures 2 and 3 provide a broad overview of the evolution of the new models we computed. 
They illustrate the time variation of energy components (thermal, kinetic and magnetic), as well 
as the pressure minimum ratio, /3mm = ip/Pb)min, on the grid at each time. Figure 2 displays 
results for the runs that were computed with Nx = 512 while Figure 3 shows the cases computed 
with Nx = 256. In addition, to provide a sense of the influence of numerical resolution for very 
weak field cases both Case 6h {Nx = 512) and Case 6m (Nx = 256) are shown together in Figure 

2. Resolution issues were discussed in detail for Cases 1 and 2 in Paper I, where similar energy 
plots were also given. We shall add a few additional relevant comments below. 

There is one distinctive detail about Figures 2 and 3 that is important to their interpretation. 
It was apparent to us by comparing animations of the important dynamical quantities that 
Cases 3 and 4 were virtually indistinguishable from one another, and likewise for Cases 6m 
and 7. Thus the important issues are somewhat easier to see if we eliminate Bz to first order 
in Figures 2 and 3. We plot a reduced energy for the cases with 6^0. The magnetic and 
kinetic energy plotted include only the planar components; i.e., E'^ = J(1/2)Bpdxdy, where 
Bp = ,Jb‘^ + B"^, and E'f, = /(l/2)p(tt^ +Uy)dxdy. To compensate, the total energy is also reduced 
as = E{j + Et + E'f., where Et is the thermal energy. Thus, we ignore the almost constant 
energy contributions from Bz and Vz. Note, however, it is the total energy, not the reduced 
total energy, that is conserved. In the definition of (drain, total magnetic pressure from all three 
components was used, since it is the entire pressure that exerts a force on the plasma. 

We can see that the reduced energy evolution in Cases 3 and 4 are almost identical in Figure 
2, and the same is true of Cases 6m and 7 in Figure 3. Thus, we observe in 2^D, at least, that the 
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transverse magnetic field component (B^) plays no significant role in nonlinear evolution of the 
instability. The rationale for some role comes from the observation that finite and enable 
circularly polarized Alfven waves in 2^0, while only linearly polarized Alfven waves are allowed 
in 2D. In other words, there are twice as many degrees of freedom in 2^D for Alfven waves to 
help disperse perturbations. However, there are two arguments to support our observation of a 
minimal role. First, the group velocity vector for Alfven waves along which physical information 


propagates is Vg = Bjyfp\ i.e., aligned with the magnetic field (see, e.g., [Landau, Lifshitz,^; 


Pitaevskii (1984)[) . With invariance of quantities along the z-direction, only the group velocity 
components projected onto the x — y plane are relevant, and they are independent of B^- The 
second point is that even though Bz and Vz are finite in 2^D, the above symmetry restricts their 
contributions to those of magnetic pressure (see equations 


|2-5|); i.e., to the total pressure 
gradient and the fast and slow wave speeds. In the present case with weak fields, even that has 
little importance. For example, the evolution of /3mm is almost the same for Cases 3 and 4, where 
stretching of the field in the plane is dominant in Case 3. For Cases 6m and 7 the values of Pmin 
are always large, because neither the planar nor the transverse fields are very strong. But at r ~ 5 
when the stretching of the planar field is greatest, the values of /3mm are still similar. So, our 
results indicate that for weak fields the extra degrees of freedom allowed in 2^D have no significant 
effect on the KH instability. These flows are essentially 2D. 


In one of the most striking findings from these simulations, we find two distinctive evolutionary 
behaviors for weak field flows, depending on the field strength. When the planar magnetic fields 
are very weak, it turns out that the evolution of the 2D MHD KH instability is qualitatively 
similar to the gasdynamic version of the instability. That is, the cat’s eye vortex continues to exist 
as long as we extend the simulation. The magnetic field does, through unsteady reconnection, 
substantially enhance energy dissipation of the vortex over that from the gasdynamical case, 
however. This can be seen by comparing Case 5 {Bp = 0) in Figure 3 with Case 6m or Case 7 
{Bp = 0.02). We will characterize such flows as having “very weak fields”, and the role of the 
field as “dissipative”. On the other hand, if the initial fields approach the critical field for linear 
stabilization within a factor of a few, the cat’s eye is disrupted and the field causes the flow to 
reorganize into a laminar form (as in Case 2). This we call simply the “weak-field” regime, and the 
field as “disruptive”. In either situation the magnetic energy peaks about the time of the initial 
reconnection instability. For the disruptive cases there is a fairly prompt return of the magnetic 
energy near to its initial value (see Figures 2 and 3). In effect, the magnetic field plays the role 
of a catalyst, storing energy temporarily and through it modifying the plasma flow. That role is 
fairly dramatic, since reconnection leads to dynamical alignment and self-organization in the flow. 
For dissipative cases, the magnetic energy declines much more slowly, and in the high resolution 
Case 6h, seems mostly to oscillate. That behavior results from the fact that the cat’s eye vortex 
continues to capture magnetic flux, amplifying it and thence dissipating it. 


The existence of two qualitatively distinct evolutionary patterns is also quite apparent in the 
histories of the thermal and kinetic energies displayed in Figures 2 & 3. The kinetic energy decay 
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in Case 6 is almost exponential, with a time constant that can be roughly estimated as ~ 10^. 
By contrast Cases 3 and 4 show before r ~ 10 a sharp drop in E'f^ along with an accompanying 
increase in E^, followed by a slow, possibly exponential evolution that is similar to Case 6. For all 
cases in Table 1 the magnetic energy is always small (although it can briefly increase by factors 
between 6 and 20 it never contributes more than a few percents to the total energy), so on the 
face of it the magnetic energy would not seem to be important. In fact, it can be crucial, as we 
have already outlined and will discuss more fully, below. Figure 3, which shows a wider range of 
projected field strengths, shows similar patterns with some variation in the abruptness of the early 
energy transition. We also conclude that as —> 0 the evolution of the 2D MHD KH instability 

will smoothly approach that for ordinary gasdynamics. 


3.1. Very Weak Fields: Dissipative 


Let us now describe in more detail the characteristics of the two weak-held regimes and 
also establish the physical boundary between them. Case 6 and Case 7 had the weakest finite 
planar field, Bpo, that we considered, and they demonstrate very-weak-field patterns. Figure 4 
summarizes flow properties for Case 6h at times r = 5 and 30. Figure 4a shows the flow vorticity 
component out of the plane, tOz = {V x u)z, while Figure 4b shows the magnetic field lines as 
they project onto the plane. Except for minor fine structure that aligns with the planar field, the 
vorticity at r = 5 is the same as for a flow with Up = 0 of Case 5. It also is fairly close to that for 
all our simulations, except Case 1 at this early time . Primarily this image illustrates the formation 
of the cat’s eye and how that concentrates vorticity. Note that the “corners” of adjacent cat’s 
eyes overlap, with a shear layer between them. It is in that shear layer where magnetic fields are 
most strongly enhanced, as seen for this time in Figure 4b. Extensions of that feature around the 
perimeter also contain concentrations of magnetic flux. These regions represent flows where gas 
is accelerated out of a stagnation point midway between vortex centers. Frozen-in fields are thus 
pulled out or “stretched” and amplified. At r = 30, the cat’s eye vortex is still largely the same, 
except for some complex, low level vorticity structures outside the main vortex. Their origins are 
made clearer by examination of the magnetic field structure at this time. Those same regions 
outside the cat’s eye also contain isolated magnetic flux islands and field reversal regions. Such 
magnetic features reveal an environment where magnetic reconnection is active. The relationship 
is that reconnective processes not only reorganize the magnetic field topology and release magnetic 
energy, but they also accelerate the local plasma, and that contributes to the local vorticity. By 
contrast to Case 2 (or as we shall see Cases 3 and 4, as well), however, magnetic stresses in Cases 
6 or 7 produce only small modifications to peripheral flows, and are far too weak to disrupt the 
cat’s eye. 


The vortex interior shows initial signs in Figure 4b of magnetic flux expulsion by r = 5, 
a well-known phenomenon {e.g., Weiss (1966) ). At this early stage, just as the vortex is fully 
formed, there is still some magnetic flux that threads through the “eye”. That does not seem to 
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be the case at the later time. In a close examination of field structures within the vortex we can 
find no evidence after about r = 10 that any magnetic flux threads the vortex. Instead, the field 
breaks into flux islands within the vortex, and those are annihilated through mergers. This is, of 
course, just what high Reynolds number (both kinetic and magnetic), nonlinear resistive MHD 
flows are expected to do in response to reconnection {e.g., Biskamp (1993) ; Biskamp (1994)|) . 
The field structure in the bottom panel of Figure 4b is qualitatively very similar to that found 
by Weiss ( Weiss (1966) ) from a classic passive field simulation in a steady vortex with magnetic 
Reynolds number Rm = 10^, but very different from that for small magnetic Reynolds number 
{Rm = 20). Our own estimates of effective magnetic Reynolds numbers in the Case 6h simulation 
give numbers in excess of 10^ ( Ryu, Jones &: Frank (1995) ; Frank et al. (199^ ), so the comparison 
is very reasonable. 


We emphasize, however, that the field behaviors seen here, and presumably by Weiss, result 
from localized, inherently time dependent reconnection, not simple flux diffusion. To demonstrate 
that we compute in the Appendix the equilibrium passive magnetic fields for a simple vortex in 
a resistive fluid. Figure 5 displays two examples of magnetic field structures predicted by this 
steady-state resistive MHD theory. The velocity structure for this model vortex is similar to 
that observed in the simulated flows for Case 6. The core, r < ro/2, has a constant vorticity. A, 
(see Figure 4a) while there is an outer flange in which the velocity decreases to zero at r = tq. 
Solutions depend only on the effective magnetic Reynolds number, q = 2ArQ/ri, within the vortex. 
This parameter roughly measures the ratio of the timescale for magnetic field diffusion to the 
rotation period of the vortex. When q is a few {q = 10 in the top panel of the figure), flux diffuses 
quickly enough to almost flatten out field lines to their external pattern except in the vortex core. 
At the center field lines are actually concentrated by rotation into a quasi-dipole pattern. That 
dipole results because reconnection into a simpler topology is not permitted. 

The steady-state pattern is very different when q is large, since field is almost frozen into the 
vortex. Above q ~ 10^ the field forms a spiral pattern in the vortex core whose pitch depends on 
q. Shear is greatest in the flange region of this vortex and flux becomes concentrated increasingly 
into the vortex perimeter. That particular behavior does resemble Case 6 in some ways. There 
are multiple field reversals in the perimeter, with that number determined by the balance between 
held diffusion and addition of new turns to the field through rotation. The bottom panel in 
Figure 5 shows the steady-state solution for q = 10^. Despite some superficial similarities, there 
is an obvious and essential difference between the truly steady-state solution in Figure 5 and the 
“quasi-steady” field in the two time dependent simulations (ours and Weiss’). For the steady 
solutions the held line that was originally through the vortex center remains and is wrapped 
into a spiral pattern within the core giving an almost constant held strength there. In the time 
dependent cases, reconnection changes the field structure inside the vortex. The same statement 
applies to Weiss’ simulation. If the effective magnetic Reynolds number within the flow is large 
enough for topologies to develop that are tearing-mode unstable, field within the vortex “breaks” 
off from the external field and is destroyed. Thus, we And good agreement with Weiss in the 
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nature of eddy flux expulsion and see that it cannot be viewed as a steady magnetic diffusion. 

The value of Bp never drops to zero inside the vortex of any of our very weak field cases, 
although it is as much as two orders of magnitude smaller than along the vortex perimeter. Action 
around the perimeter episodically injects new magnetic flux into the vortex interior, so this process 
continues as long as we have followed the flow. Examination of the magnetic energy, E'^, for Case 

6 in Figure 2 shows related, episodic peaks on rough intervals, At ~ 2.5 — 3. That corresponds 
to about half a turn over-time for the vortex, and represents the interval on which fields along the 
vortex perimeter are stretched until they become “folded”, so as to make them subject to tearing 
mode reconnection. So, whatever flux is caught in this flow becomes amplified, increasing the 
magnetic energy, before reconnection rearranges the field lines. Those lines outside the vortex are 
relaxed towards the initial field configuration, often through multiple reconnection events, while 
those inside are isolated into closed islands. Note that the total magnetic flux threading our box 
is exactly conserved and exactly zero in these computations. Thus any field line entering on the 
left boundary must at all times extend continuously to the right boundary or exit again on the 
left. Closed flux islands can exist, but nothing prevents them from being destroyed, since they 
contribute no net flux. 

Reconnection, and most obviously flux-island destruction, is irreversible, so that it must 
be accompanied by entropy generation. That outcome is very apparent in Figures 2 and 3. 
Recall that to first order total thermal energy changes in these, closed-system simulations reflect 
nonadiabatic processes. Compare first the evolution of thermal energy in Cases 6m, 7 and 5. We 
see that after some viscous dissipation necessary to form the cat’s eye, the Case 5 with BpQ = 0 
has almost constant total thermal energy to the end of the simulation. By contrast Cases 6m and 

7 show a steady rise in Et, with small amplitude oscillations associated with major reconnection 
events as described above. Thus, the main dynamical impact of the field is enhanced dissipation. 
The total magnetic energy remains small in all cases, but for the dissipative cases that reflects 
a near balance between the rate at which kinetic energy is being transferred to magnetic energy 
and the rate at which magnetic energy is being dissipated. For example, from initial field energy 
evolution in Figure 3 we can crudely estimate for Case 6m that magnetic energy is generated at a 
rate d E'^^/dt ~ 3 x 10“^, which is very close to the mean slope of the thermal energy curve, Et- 

Looking next at the Case 6m,h Et plots in Figure 2, we see that energy dissipation is greater 
in the higher resolution Case 6h. That results despite a smaller numerical resistivity for the higher 
resolution simulation. This is because in Case 6h a larger amount of flux is caught in the vortex, 
so more reconnection sites develop. We estimate that ~ 10% of the total magnetic flux is attached 
to the vortex at the end of the calculation (r = 50) in Case 6h, while only ~ 3% is attached at the 
end of the calculation (r = 30) in Case 6m. As the effective magnetic Reynolds number increases 
there is a tendency for decreased resistivity to be countered by an increase in the number of 
reconnection sites, as our earlier discussion of reconnection theory would suggest. If we compare 
simple estimates for rates of magnetic energy generation for Cases 6m and 6h to the rates of 
thermal energy increase, we see that they are consistent, as we found earlier for Case 6m, alone. 
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The increased dissipation for the higher resolution simulation is relatively modest, however. Some 
studies of resistive MHD turbulence suggest, in fact, that at very large Reynolds numbers total 
reconnection rate in a complex flow will be insensitive to the value of the resistivity (Biskamp 
(1993); Biskamp (1994)). If confirmed more generally, that could provide a practical measure of 


convergence in studies of the present kind. 


3.2. Weak Fields: Disruptive 

Looking again at thermal energy evolution in Figures 2 and 3 we can see a clear behavioral 
transformation in the sequence: Case 5 —> 6(7)—> 8 —> 9 —> 3(4). This grouping is arranged in order 
of increased BpQ. There is a sharp rise through the sequence in the amount of dissipation associated 
with the initial formation of the cat’s eye vortex, as well as a more modest increased slope to the 
subsequent long-term dissipation rate. The physical character of this transformation is apparent 
if one compares Figure 4b with Figure 6. The latter shows for Case 3 the magnetic pressure 
distribution at three times, r = 5,10 and 30. Behaviors for this simulation are qualitatively similar 
to Case 2, as discussed in detail in Paper I and outlined in §2.3. Here we note that, while the field 
appears wrapped around the vortex at r = 5, it has a laminar appearance at r = 30. The velocity 
field undergoes a similar transition; i.e., the vortex is completely disrupted. The magnetic field at 
the intermediate time, r = 10, shows aspects of both the other times. Curiously, however, the flow 
in the dominant vortex pattern there has the opposite vorticity to the original flow. That feature 
is short-lived. We can understand this flow transition from rotational to laminar by examining 
Figures 7, 8 and 9, also relating to Case 3. In Figure 7 we show the log of /3 = p/pb with magnetic 
field lines overlaid at r = 5 . The minimum /3 = 0.55 in the strong flux tube connecting vortices, 
but the magnetic field is dynamically significant most of the way around the perimeter of the 
vortex. The maximum f3 ~ 10® in small regions where magnetic reconnection has begun and the 
field strength has decreased to very small values. 

To understand how the magnetic field disrupts the vortex consider the forces involved. The 
centripetal force associated with motion around the cat’s eye is pu^/R hpui ~ since we 
observe that ~ i? ~ 1 for the vortex. Indeed we also confirm that when BpQ = 0, the 

pressure gradient force within the vortex, —Vp ~ 0.5 and is directed towards the vortex center, so 
that it supplies the necessary force. On the other hand in Figure 8 we display the magnetic field 
lines along with vectors representing the magnetic tension force; i.e., 

T=(B-V)B, (3-1) 

at r = 7.5 for Case 3. It is apparent that magnetic tension forces are concentrated where the 
field has been pulled into loops by the flow, and that they are directed towards the center of the 
vortex. The peak value of |T| = 1.85. At r = 7.5 in Case 3, however, the total pressure gradient 
is small and actually of the wrong sign to effect significantly the motion of the plasma. The 
total force vector field is very similar to the tension force field shown in Figure 8. At the same 
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time “X-points” in the field topology show that reconnection is underway that will isolate the 
associated field lines. Subsequently, flux islands are formed along the “axis” of the cat’s eye, and 
the magnetic tension pulls the plasma frozen into those loops towards the original center of the 
vortex, disrupting its rotation. Field line segments reconnecting outside the vortex core will tend 
to relax towards the original field topology. 

These observations allow us to estimate simply what minimum initial field, Bp^ should lead to 
vortex disruption. Since the early evolution of a weak field is self-similar, we can use the behavior 
from Case 3 to estimate |T| ~ 1.85 (i3po/0.14)^ at the time of the first major reconnection in all 
weak-field cases. Our earlier discussion requires |r| > \UqP = ^ for disruption, leading to the 
constraint Bpo > 0.05, or alternately, Bpo > O.lBc, where Be is the critical field for stabilization 
of original instability. Indeed, as Figure 3 demonstrates, the transformation between dissipative 
and disruptive evolution occurs for conditions between those of Cases 8 and 9; i.e., for initial field 
values between 0.04 and 0.07. 

Figure 9 illustrates why the vortex disruption process is also highly dissipative. The top 
panel shows at r = 7.5 in Case 3 the gas entropy distribution, while the lower panel displays 
the electric current density, \j\ with the field lines overlaid. From this we can see that excess 
entropy is concentrated into regions where reconnection is currently active (highlighted by Ijj) or 
recently active. Ohmic heating (oc is partly responsible for the irreversible energy exchange. 
The remainder should be viscous dissipation of small-scale, disordered motions. As we discussed 
in Paper I, the final laminar flow that results in this class of flow includes a central sheet of hot 
gas containing most of the excess entropy produced through the self-organization of the flow. 


4. Summary and Conclusion 

We have carried out a series of high resolution MHD simulations of Kelvin-Helmholtz unstable 
flows in 2^ dimensions. All of these simulations involve magnetic fields initially too weak to 
stabilize the flows in the linear regime; i.e., BpQ < Be. Thus, since simulations are performed on 
a periodic space, flows all begin formation of a single “cat’s eye” vortex. If the field lying in the 
computational plane is absent or “very weak” the cat’s eye structure becomes a persistent, stable 
feature that represents a “quasi-steady” equilibrium. When there are “very weak” magnetic fields 
in the plane they become wrapped into the vortex and amplified by stretching. However, within 
a single turn of the vortex they are subject to tearing mode instabilities leading to magnetic 
reconnection. That reconnection isolates some magnetic flux within the vortex, which is eventually 
annihilated. This is the process through which flux is effectively expelled from a vortex. As long 
as the vortex persists this process will repeat. Since reconnection is irreversible, this process is 
also dissipative and leads to an increase over viscous effects in conversion from kinetic to thermal 
energy. We find in this regime that as the initial magnetic field within the computational plane is 
increased the dissipation rate increases in a similar manner. Likewise, as we use a finer numerical 
grid, thus reducing the eff'eetive numerical resistivity and viscosity, the dissipation rate inereases, 
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reflecting the increased ability of our code to capture small-scale reconnection events. This trend 
is backwards from what one would expect if simple magnetic diffusion were primarily responsible 
for the reconnection. It suggests, perhaps that if we had been able to extend these calculations to 
even higher resolution the energy dissipation rate might have converged to a value independent of 
the effective resistivity, just as some studies of resistive MHD turbulence find. The reconnection 
and expulsion of flux within vortices in our simulations are similar to those in a classic study by 
Weiss of vortex flux expulsion in large-magnetic-Reynolds-number flows. 

If the initial magnetic field is strong enough that within a single turn of the vortex it is 
amplified around the vortex perimeter to “dynamical” strength (5^ ~ pU"^)-, then the reconnection 
described in the previous paragraph releases stresses that are capable of disrupting the vortex 
entirely. This can happen in a single event or, if the field is only marginally strong enough 
{BpQ through a succession of dynamical realignment events. In either case the net result 

is a laminar, but marginally stable flow, in which the original shear layer is greatly broadened. 
Thus, as we discussed fully in Paper I, such fields can have a remarkable stabilizing influence. 
This is despite the fact that their total energy content is a minor fraction of the total, so that they 
are nominally too weak to be important, according to the usual criteria. 


We considered cases in which the magnetic field was entirely within the flow plane and others 
in which the field was oblique to that plane, in order to examine the role in nonlinear flows of 
the component out of the plane. For the 2^ D flows we have studied, only the field components 
in the flow plane have any dynamical significance. In fully 3 D flows, however, we expect further 
evolution of the “quasi-steady relaxed states” of both very weak field (or dissipative) cases and 
weak field (disruptive) cases. The cat’s eye vortex of very weak field cases is subject to a 3 D 
instability known as the elliptical instability ( Pierrehumbert (1986) ; Bayly (1986) ) unless the flow 
lines around the vortex follow perfect circles. The planar shear flow of weak field cases is stable 
against linear perturbations but unstable to 3D finite-amplitude perturbations (Bayly, Orszag & 
Herbert (1988) ). Thus, it will be important to extend the present study to the fully 3 D regime, 
and we are preparing to do that. 


This work by TWJ, JBG, and AF was supported in part by the NSF through grants 
AST-9318959 and INT-9511654, through NASA grants NAGW-2548 and NAG5-5055 and by the 
University of Minnesota Supercomputer Institute. The work by DR was supported in part by 
Seoam Scholarship Foundation. We are grateful to B.I. Jun for stimulating and helpful discussions 
about these results. 


A. Diffusive Flux Expulsion from a Steady Vortex 

In order to show that the flux expulsion from vortex shown in Figure 4b requires localized, 
unsteady reconnection, not just diffusion, here we study the steady state passive field solutions 
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of resistive MHD. Under 2D and 2^D symmetries, the planar magnetic field can be written as 
B = z X S/tp where xp is the magnetic flux function. Then, the induction equation in (2-5) with 
constant resistivity becomes 

f)-,h 

(Al) 


^ + u ■ Vxp = ■qV'^'ip. 
ot 


We consider the evolution of an initially uniform magnetic field, B = BqX ov xp = 
within a vortex. We approximate that the vortex has azimuthal velocity 


2Ar, for 0 < r < ro/2 

v,j)= \ 2A{ro - r), for ro/2 < r < ro 
0, for r > ro 


(A2) 


and zero radial velocity. A equals the (constant) vorticity within the vortex core. This is roughly 
the velocity field of the vortex in Figure 4. Then, the evolution of the magnetic field is described 
by equation ([Al|) with the boundary conditions xp{ro) = BoVoC^^ and dxp{ro)/dr = Boe^'^ at r = Tq. 
In a steady state, we can set xp = F(r)e*"''^ since the coefficients of the equation do not depend 
explicitly on (p. Note that with the given initial magnetic field and boundary conditions, only the 
solution with n = 1 is allowed. Then, the equation for F is given as 


dF I dF 1 . Q p _ 

r dr ^ r^ 

where 

{ q, for 0 < r < ro/2 

q{ro/r — 1), for ro/2 < r < ro 
0, for r > ro 

with q = 2Ar/^/r]. Here, q is the magnetic Reynolds number within a factor of two. 


(A3) 


(A4) 


We solve equation (A3) for F numerically. In Figure 5, we plot the resulting magnetic field 
lines for the case with g = 10 (highly diffusive case) and for the case with q = 10^ (quasi-adiabatic 
case). From dissipation tests we estimate the magnetic Reynolds number of the vortex in our 
time-dependent simulations to be larger than 10^ (see Ryu, Jones Sz Frank (I995)| ; Paper I). 
Indeed the fields in Figure 4b resemble those in a vortex simulation with magnetic Reynolds 
number 10^ reported in Weiss (I966)| . The steady-state solution we just described allows only for 
diffusion of magnetic flux, since no change in the field topology is permitted. It allows us to see 
that such a steady-state, diffusive description does not account for flux expulsion from the vortex. 
If that were the case the magnetic field lines in Figure 4b should have a structure like those in 
the bottom of Figure 5. Actually, their topologies are fundamentally different in the sense that 
field in Figure 5 threads completely through the vortex, while it does not in Figure 4b (or any 
similar figure showing additional field lines for this simulation). The reason is that the field in the 
time dependent simulation is subjected to reconnective instabilities that isolate and then destroy 
magnetic flux in the interior of the vortex. 
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Table 1. Summary of MHD KH Simulations 


Case 

Bo = Ca/Cs'" 

e 

Bpo = Bo cos 0 

Ma'° 


f ^ 

End Time° 

iV/ 

1 “ 

0.4 


0.4 

2.5 

7.5 

3.79 

20t 

512“ 

2 “ 

0.2 


0.2 

5.0 

30 

1.86 

20 r 

512“ 

3 

0.2 

45° 

0.14 

5.0 

30 

1.71 

30r 

512 

4 

0.14 

0 ° 

0.14 

7.07 

60 

1.71 

30r 

512 

5 

0.2 

0 

0 

0^ 

0.0 

5.0 

30 

1.59 

20 r 

256 

6 h 

0.2 

0 

00 

0.02 

5.0 

30 

1.59 

50r 

512 

6 m 

0.2 

85° 

0.02 

5.0 

30 

1.59 

30r 

256 

7 

0.02 

0 ° 

0.02 

50 

3000 

1.59 

30r 

256 

8 

0.04 

0 ° 

0.04 

25 

750 

1.67 

30r 

256 

9 

0.07 

0 ° 

0.07 

14.3 

245 

1.67 

20 r 

256 


^All models have used 7 = 5/3, M = Uq/cs = 1, Cs = 1, L = 2.51, a — L/25. 

’^The Alfven speed, Ca = B /Alfven Mach number, and /3o = Pg/Pb = ( 2 / 7 )(Ma/M)^ here refer to the total 
initial magnetic field strength, not just that projected onto the plane of the flow. 

'^The growth time is an approximation to the inverse linear growth rate; namely tg = l/F; t = t/tg. 

‘^Computations were carried out on a square grid of the size indicated, with Ny = N^- Periodic boundaries were 
assumed for x and reflecting boundaries were assumed for y. 

“Cases 1 and 2 were presented and discussed in Paper I. They are cited here for reference. Each was computed 
with 2 or more numerical resolutions, with the largest listed here. 
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FIGURE CAPTIONS 

Fig. 1 Cartoon illustrating the computational setup for these simulations. There is a central 

shear layer separating uniform flows. Except for the initial perturbation the flow conditions 
are otherwise uniform. The magnetic field projects onto the computational plane at an 
angle, 9. 

Fig. 2 Time evolution of the high resolution simulations (Cases 3, 4, 6 h), plus the medium 
resolution simulation Case 6 m. Shown are the normalized total thermal, kinetic and 
magnetic energies, as well as the minimum value of the plasma /3 parameter at each time. 
To emphasize the minimal importance of the transverse magnetic field component, Bz, 
“reduced” energies are shown that exclude Bz- 

Fig. 3 Same as Figure 2 , except for medium resolution simulations (Cases 5, 6 m 7, 8 , 9). 

Fig. 4a Greyscale images of the Case 6 h vorticity component normal to the computational plane 
at r = 5 (top) and r = 30 (bottom). To facilitate visualization of structures, the periodic 
space has been repeated once. The vorticity is everywhere positive. High values take high 
tones. 

Fig. 4b Case 6 h magnetic field lines projected onto the computational plane plane at r = 5 (top) 
and r = 30 (bottom). 

Fig. 5 Magnetic field lines from an analytical model to study diffusive flux expulsion from a 
steady eddy with (7 = 10 (top) and q = 1000 (bottom). See the text for the definition of q. 

Fig. 6 Grayscale snapshots of the magnetic pressure distributions for Case 3 at r = 5, 10 and 30, 
showing the stages in disruption of the cat’s eye vortex by the magnetic held. High values 
take high tones. 

Fig. 7 Inverted grayscale image of log/3 = p/pb for Case 3 at r = 5. Projected magnetic held 
lines are overlaid. The minimum /3 ~ 0.55 in the strong hux tube between vortices, while the 
maximum /3 ~ 10® near the center of the cat’s eye. 

Fig. 8 The magnetic held structure for Case 3 at r = 7.5, just as substantial magnetic 

reconnection is underway. Projected held lines are shown, along with arrows that represent 
magnetic tension forces. The maximum magnetic tension force is 1.74 

Fig. 9 Inverted grayscale image of the gas entropy (top) and electrical current (bottom) 

distributions for the Case 3 at r = 7.5. Magnetic held lines are also laid on top of the current 
distribution to emphasize the relationships. High values take low tones. 
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